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Abstract 



^ , We present a new high-quality nucleon-nucleon potential with explicit charge 

h : 

dependence and charge asymmetry, which we designate Argonne vis- The 
model has a charge-independent part with fourteen operator components that 
is an updated version of the Argonne t»i4 potential. Three additional charge- 
dependent and one charge-asymmetric operators are added, along with a com- 
plete electromagnetic interaction. The potential has been fit directly to the 
Nijmegen pp and np scattering data base, low-energy nn scattering parame- 
ters, and deuteron binding energy. With 40 adjustable parameters it gives a 
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X 2 per datum of 1.09 for 4301 pp and np data in the range 0-350 MeV. 
13.75.Cs, 12.40.Qq, 21.30.+y 
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I. INTRODUCTION 



Traditionally, nucleon-nucleon (AW) potentials are constructed by fitting np data for 
T = states and either np or pp data for T = 1 states. Examples of potentials fit to np data 
in all states are the Argonne vu [jl], Urbana Vu @, and most of the Bonn potentials 0,f|. 
In contrast, the Reid ||, Nijmegen ||, and Paris M potentials were fit to pp data for T = 1 
channels. Unfortunately, potential models which have been fit only to the np data often give 
a poor description of the pp data [Q , even after applying the necessary corrections for the 
Coulomb interaction. By the same token, potentials fit to pp data in T = 1 states give only a 
mediocre description of np data. Fundamentally, this problem is due to charge-independence 
breaking in the strong interaction. 

In the present work we construct an updated version of the Argonne potential that fits 
both pp and np data, as well as low-energy nn scattering parameters and deuteron properties. 
The strong interaction potential is written in an operator format that depends on the values 
of S, T, and T z of the AW pair. We then project the potential into a charge-independent 
(CI) part that has fourteen operator components (as in the older Argonne Vu model) and 
a charge-independence breaking (CIB) part that has three charge-dependent (CD) and one 
charge-asymmetric (CA) operators. We also include a complete electromagnetic potential, 
containing Coulomb, Darwin-Foldy, vacuum polarization, and magnetic moment terms with 
finite-size effects. We designate the new model Argonne v\g. 

In a number of applications it is important for an AW potential to reproduce correct np 
and pp scattering parameters. For example, in thermal neutron radiative capture on the 
proton, p(n, 7)d, it is crucial to have the correct singlet np scattering length in the initial 
state to get the cross section. However, in low-energy proton weak capture, p(p,e + b , e )d, it 
is equally important that the correct pp scattering length be provided by the interaction. 
Clearly, a complete potential model should meet both requirements. 

Another important application is in the formulation of three-nucleon (NNN) potentials. 
In general, nuclei are underbound using only AW potentials fit to the scattering data. Non- 
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trivial many-nucleon interactions are expected to make up a portion of the missing binding 
energy. Phenomenologically we may choose to construct a many-body Hamiltonian, such as 

# = E^ v *+£% + £ V ijk , (l) 

i ZIVI i i<j i<j<k 

and constrain the strength parameters of the NNN potential by requiring that H gives the 
correct trinucleon binding energy. Similar considerations apply if we choose a relativistic 
formulation. Clearly, such constraints are ambiguous or even meaningless if the NN potential 
used in the calculations does not adequately describe the two-nucleon data. For 3 He ( 3 H), 
in which the NN interaction underbinds by ~1 MeV, there are two np pairs and one pp 
(nn) pair. To a good approximation, the two np pairs will be in the S = 1, T = state 
75% of the time, and in the S = 0, T = 1 state 25% of the time, while the pp (nn) pair 
will be pure S = 0, T = 1. If the chosen NN potential fits only the more repulsive pp (nn) 
data in the T = 1 state, we would get a smaller NN contribution to the binding energy and 
thus overestimate the NNN potential strength required. By the same token, a model fit to 
np data in the T = 1 state would be too attractive and we would underestimate the NNN 
potential. The difference can be as much as 0.4 MeV, leading to variations in the NNN 
potential strength of order ±20%. This would have significant effects in larger many-body 
systems. 

Because we include a complete electromagnetic potential and fit low-energy nn scattering 
data, the present model also can be used to study charge- symmetry breaking, as in the 3 H - 
3 He mass difference ||, or more generally the Nolen-Schiffer anomaly |10 |. The electromag- 



netic potential is in principle well-known and is the longest-range part of the interaction. 
Potential models commonly fit the deuteron energy to better than 1 keV accuracy. Since we 
find that the electromagnetic terms give a non- negligible 18 keV repulsion in the deuteron 
and moderate shifts in the np and nn scattering lengths, we deem it desirable to include 
these terms explicitly. 

The major goal of the present work is to construct a nonrelativistic potential that can be 
used easily in nuclear many-body calculations and that accurately fits both pp and np data. 



We adopt the local operator structure of the older Argonne vu and Urbana Vu potentials, 
which have been used extensively in calculations of finite nuclei, nuclear matter, and neutron 
stars pH|JT^JT3[1 . The assumption of an underlying operator structure relates all partial waves 
in a simple manner, without imposing a one-boson-exchange (OBE) form which might be 
too restrictive at short distances. Recently, the Nijmegen group has shown |TJJ that it is 
feasible to construct potential models which fit the NN data with the almost perfect x 2 P er 
datum of 1. However, these models differ in each partial wave and thus implicitly introduce 
nonlocalities from one partial wave to the next that may be difficult to characterize and treat 
accurately in many-body calculations. When they limit the potential to an OBE form, which 
has a local operator structure (save for a nonlocal part in the central potential) describing all 
partial waves simultaneously, the x 2 P er datum increases to 1.87, albeit with a much smaller 
number of parameters. The present model is a compromise between these two approaches, 
adopting a phenomenological form (unrestricted by an OBE picture) at short distances, but 
maintaining a local operator structure. The potential was directly fit to the Nijmegen NN 
scattering data base [|T^,|T^], which contains 1787 pp and 2514 np data in the range 0-350 
MeV, and has an excellent x 2 P er datum of 1.09. 

In Sec. [TT] we present the analytical form of the potential in the various spin and isospin 
states. Special attention is given to the electromagnetic part of the interaction. The free 
parameters are fit to the NN scattering data and deuteron binding energy in Sec. [TTJ, where 



we also present the phase shifts. Section [TV] discusses the projection of the potential into 
operator format. Static deuteron properties and electromagnetic form factors, with rela- 
tivisitic and exchange current contributions, are presented in Sec. [V[ Conclusions and an 



outlook are given in Sec. VI. 
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II. FORM OF THE POTENTIAL 
IN S, T, T z STATES 



The NN potential is written as a sum of an electromagnetic (EM) part, a one-pion- 
exchange (OPE) part, and an intermediate- and short-range phenomenological part: 



v(NN) = v (NN) + v n (NN) + v (NN) 



(2) 



The EM interaction is the same as that used in the Nijmegen partial-wave analysis, with 
the addition of short-range terms and finite-size effects |IT||T8| , |l9"f . (Values for the masses 
and other physical constants used in the following formulae are given in Table |I[) For pp 
scattering we include one- and two-photon Coulomb terms, the Darwin- Foldy term, vacuum 
polarization, and the magnetic moment interaction, each with an appropriate form factor: 



v EM (pp) = V C i(pp) + V C2 + V DF + V VP + V MM (pp) ■ 



(3) 



Here 
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The Coulomb interaction includes an energy dependence through the a' = 2ka / (MpV^) p0| . 
which is significantly different from a at even moderate energies (~ 20% difference at Tj a b = 
250 MeV). The vacuum polarization and two-photon Coulomb interaction are important for 
fitting the high-precision low-energy scattering data. The Fq, F§, F t , and Fi s are short- 
range functions that represent the finite size of the nucleon charge distributions. They have 
been obtained under the assumption that the nucleon form factors are well represented by 
a dipole form 



G p E = ^ = ^ = G D =(l + pi' 2 , (9) 



where b = 4.27 fm 1 . The functions are given by 

11 3,1 



F c (r) = l-(l + -x + -x 2 + -x 3 ) e -*, 

F t (r) = 1 - (l + x + l -x 2 + l -x z + ^x 4 + ^x 5 ) e" x , 

F is (r) = 1 - (l + x + \x 2 + ^x 3 + ^x 4 ) e~\ (10) 



with x = br. The derivation of _Fc is given in ref. ||21|| , while the others are related by 
F 5 = -V 2 (F c /r), F t = {F c /r)" - (F c /r)'/r, and F ls = (F c /r)'/r. In the limit of point 
nucleons, Fq = F t = F[ s = 1 and Fg = 47r<5 3 (r). These form factors are illustrated in 
Fig. |l| The use of Fq in Vyp is an approximate method of removing the 1/r singularity 
(the logarithmic singularity remains) which is justified by its short-range and the overall 
smallness of the term. Similarly, the use of Fq in Vc% is an approximate method of removing 
the 1/r 2 singularity. We note that because we use the Sachs nucleon form factors, there are 



no additional magnetic Darwin- Foldy terms [22 



For the np system we include a Coulomb term attributable to the neutron charge distri- 
bution in addition to the interaction between magnetic moments, 

v EM (np) = V cl (np) + V MM (np) . (11) 

Here 

Voxinp) = a(3 n ^- , (12) 



where the function F np is obtained assuming the neutron electric form factor |22|j 

2\ -3 
1_ 

b 2 



G E = f3 n q 2 + • (13) 



Here (3 n = [dG E /dq 2 ] q=0 = 0.0189 fm 2 , the experimentally measured slope |2"Bfl . We have 
checked this form factor in a self-consistent calculation of the deuteron structure function 



A(q 2 ) used to extract G% [24[] and find it gives a fairly good fit to the data. This simple 
form leads to 



F np (r) = b 2 (l5x + 15x 2 + Qx 3 + 



x 



384 



(14) 



The F np is also shown in Fig. [T| The magnetic moment interaction is given by 

Vmm{up) ■■ 



a 



AM n M p 



-F s {r)(T i -(T j + —^S ij 



li n — — (L-S + L-AJ 



2M n M r 



(15) 



where M r is the nucleon reduced mass. The term proportional to A 



o\- is a 



"class IV" charge-asymmetric force |[25|| , which mixes spin-singlet and spin-triplet states. Its 
contribution is very small, and we only include it when we construct the magnetic moment 
scattering amplitude [|19[ . 



Finally, for nn scattering, we neglect the Coulomb interaction between the neutron form 
factors, so there is only a magnetic moment term 



v EM (nn) = V MM (nn) 



a 



AM 2 



-Fs{r)<T V (T j + —^-S k 



(16) 



The charge-dependent structure of the OPE potential is the same as that used in the 
Nijmegen partial-wave analysis and reads 



v T (nn) = f 2 nV7T (m n o) , 



\T+1 Q f 2 



2/> 7r (m 7r ±) 



(17) 



where T is the isospin and 



(18) 



(Strictly speaking, the neutron-proton mass difference gives rise to an OPE "class IV" 
force as well, which again we only explicitly include when we construct the OPE scattering 
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amplitude [0.) Here Y^{r) and T M (r) are the usual Yukawa and tensor functions with the 
exponential cutoff of the Urbana and Argonne V\& models 



W = — (i 

//r V 



e 



/ir 



w^V^)^ 1 -^') 2 . (19) 

where /i = mc/Ti. The scaling mass m s , introduced in Eq. ( |18D to make the coupling 
constant dimensionless, is taken to be the charged-pion mass, m^±. The Nijmegen partial- 
wave analysis of AW scattering data below 350 MeV finds very little difference between the 



coupling constants P6| , so we choose them to be charge-independent, i.e., f pp = —f n 
f c = /, with the recommended value f 2 = 0.075. Thus all charge dependence in Eqs. ([H 
is due simply to the difference in the charged- and neutral-pion masses. 

The remaining intermediate- and short-range phenomenological part of the potential is 
expressed, as in the Argonne Vu model, as a sum of central, L 2 , tensor, spin-orbit, and 
quadratic spin-orbit terms (abbreviated as c, 12, t, Is, ls2, respectively) in different S, T, and 
T z states: 

Vst(NN) = v c STjNN (r) + v l2 T)NN (r)L 2 + v^^r) S 12 + v l s s TjNN {r)L-S + v l s s2 iNN (r)(L-S) 2 . 

(20) 

Each of these terms is given the general form 

v ST,NN( r ) = ^ST,NN^( r ) + ^ST,NN + A* r QsT,NN + (/ ir ) 2 -^5T,AW W(r), (21) 

where /i = ^(m^o + 2m 7r ±)c//i is the average of the pion masses and T M (r) is given by 



Eq. flT9"|). Thus the T 2 (r) term has the range of a two-pion-exchange force. The W(r) is a 
Woods-Saxon function which provides the short-range core: 

W(r) = [l + e (r - ro)/a ] _1 . (22) 

The four sets of constants I l s TtNN , Pst,nni Qst,nni an d Rst,nn are parameters to be fit to 
data. However, we also impose a regularization condition at the origin which reduces the 
number of free parameters by one for each v l STNN . We require that 



V S T,NN( r = 0) = 



UU ST,N N 

dr 



\r=0 



. 



(23) 



Since the tensor part of the OPE potential already vanishes at r = 0, the first condition is 
satisfied by setting Pst,nn = 0- The second condition is equivalent to fixing, for i ^ t, 

1 



Q 



ST,NN 



fxW(0) 



P. 



ST,NN 



dW dvl T 



dr 



dr 



(24) 



r=0 



where we only have to evaluate the derivative of the spin-spin part of the OPE potential. 



III. DATA FITTING 

An initial survey of possible potential forms was made by fitting to the \ 2 hypersurface 
of the Nijmegen partial- wave analysis of pp and np data ||16|| . These studies helped select the 
final form of the potential (~10 variations were tried) and the values of the function shape 
parameters c, r , and a. Eventually, the cutoff parameter in the OPE functions Y^(r) and 
T M (r) was set at c = 2.1 fm -2 , while the parameters in the short-range Woods-Saxon W(r) 
were set at r = 0.5 fm and a = 0.2 fm. This value of c is slightly different from the 2.0 
fm -2 used in the Urbana and Argonne Vu models, while ro and a are the same. Attempts 
to make a softer-core model led to a poorer fit. Sensitivity to the OPE coupling constant 
was also checked before the recommended value [p6| , f 2 = 0.075, was adopted as optimal. 



Once these four parameters were set, a preliminary fit of the remaining parameters 
IstnniPstnn-iQstnn-i an d Rst nn t° the phase shifts was made. The final values were 
obtained by a direct fit to the Nijmegen pp and np scattering data base and the deuteron 
binding energy. We use nonrelativistic kinematics, i.e., the deuteron binding energy is 
taken as Ed = n 2 /2M r . In practice, we found no benefit to including an R l ST NN in spin- 
singlet states, so these values were set to zero. Also, we found no indication of a need 
for charge dependence in the phenomenological part of spin-triplet states. In the final fit 
there are 40 nonzero intermediate- and short-range parameters. In addition, we fit the nn 



scattering length and effective range as determined by d(ir } ^)nn experiments |27|, which 
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seem most consistent with charge-symmetry breaking in the trinucleon. This is done by 
slightly modifying the short-range S = 0, T = 1 pp potential, resulting in a small difference 
between Poi tPP and Pqi nn- The same difference was carried over to the Pf x nn parameter, as 
discussed below in Sec. [IV[ All these parameter values are given in Table [TT]. 

The Nijmegen NN scattering data base [|I5|,[nj includes 1787 pp data (1656 observables 
and 131 normalization data) and 2514 np data (2366 observables and 148 normalization data) 
in the ranee 0-350 MeV. The total x 2 of the potential is 4675, divided into x 2 (pp) = 1962 
and x 2 ( n P) — 2713. A detailed breakdown of the x 2 , analogous to Ref. |§ for the pp data, 



is given in Table III. We also show the breakdown for the Nijmegen combined partial- wave 



analysis fll~6 |. The difference between these starts to increase beyond ~150 MeV. We should 



mention that there are a number of groups of np total cross section data which extend over 



a wide energy range. So in order to present the results in the form of Table |T|, we had to 
split each of these groups into a number of subgroups, each contributing in its appropriate 
energy bin. Whenever one of these groups has a normalization error, we choose to apply 
this same normalization for each of its subgroups. As a consequence, the number of np data 
in Table |TTT| is increased by 12, while the total x 2 { n P) is lowered by 28. The reason for 
this reduction in x 2 is that these 12 extra normalizations are optimized for each subgroup 
separately. 

As an independent test, we have also checked our results with the scattering analysis 
interactive dial-in (SAID) program, of the Virginia Polytechnic Institute and State Univer- 
sity (VPI&SU) group f28fl . We provided SAID with phase shifts calculated at the seventeen 
energies 1, 5, 10, 25 (25) 350 MeV. The SAID program then uses an interpolation routine 
to provide the phase shifts at all energies at which the experimental data were measured. 
For the S waves below 25 MeV, this interpolation deviates slightly from the actual values of 
the potential. Moreover, the treatment of the electromagnetic contributions is less sophis- 
ticated in the VPI&SU analysis, which leads to a large discrepancy for laboratory energies 
below 2 MeV. When we then compare with the data between 2 and 350 MeV, we obtain a 
X 2 {pp) — 2107 for 1644 pp data and a x 2 { n P) = 4157 for 3020 np data, all from their data 
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set NN943. The total \ 2 per datum by this comparison is still a very respectable 1.34. 

To demonstrate once more the importance of fitting to both pp and np data, we used the 
np version of the Argonne V\s potential, included the proper electromagnetic interaction, 
and confronted it with the pp data. In this way we partially simulate the comparison with 
the pp data of a potential model fit only to the np data. Of course, the analogue is not 
perfect, because in the Argonne vi$ potential the 12, t, Is, ls2 parameters in the triplet T = 1 
partial waves are fit to both pp and np data. Still, this modified np potential gives a \ 2 P er 
datum of 4.4 for the pp data between 2 and 350 MeV, which is much worse than the \ 2 P er 
datum of 1.1 for the actual pp Argonne v\$ potential on the same energy interval. Similarly, 
we can replace the T = 1 np part of the Argonne t> 18 potential by the T = 1 pp part and 
confront it with the np data. This modified potential then gives a x 2 P er datum of 1.8 on 
the np data between 2 and 350 MeV, rather than 1.1 for the actual np Argonne t>ig. 

The L = phase shifts are calculated using the potentials discussed in Sec. [T| i.e., 
including the complete electromagnetic interaction, and by matching to electromagnetic 
wave functions (S§m+n m ^ ne notation of Ref. |2!|). For L ^ we use the fact that &em+n 
can be reasonably approximated pP|JT^[ ] by only including the Coulomb interaction with a' 



(in case of pp scattering) or no electromagnetic interaction at all (in the case of np or nn 



scattering). The resulting phase shifts for partial waves with J < 3 are shown in Table [TV 



for pp scattering, in Table |V] for nn, in Table |V| for np in T = 1 states, and in Table |VII| for 



np in T = states. Note that the non-S" pp phase shifts in Table [iy| are calculated including 
the form factor Fc{r) in the Coulomb potential Vc\{pp)- 

In addition, we show figures of some of the more interesting phases and compare to the 
Nijmegen multienergy partial- wave analysis ||16|| , the single-energy analysis from SAID JJS 



and recent single-energy analyses by Bugg and Bryan [SDj, and by Henneck In Fig. |2| 
the pp, nn, and np 1 Sq phases of Argonne v i§ are shown, and seen to be in good agreement 
with the various analyses. The charge dependence is clearly evident; a discussion of the 
relative size of various contributions to charge-independence breaking is given below in 



Sec. [TV]. The Pq phases are shown in Fig. y; the Po channel displays the second greatest 
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amount of CIB after the 1 So channel. Again there seems to be reasonable agreement with 
the various partial-wave analyses. The t\ mixing parameter, shown in Fig. |], is both the 
most difficult to determine in single-energy analyses, as indicated by the range of values and 
size of error bars, and one of the most important because of its relation to the strength of 
the tensor interaction. The Argonne v is value tracks the Nijmegen multienergy analysis up 
to Ti a b = 100 MeV before deviating slightly on the high side. However, the differences with 
the Nijmegen multienergy analysis are still within two standard deviations. Finally, the x Pi 
phase shift, which is intimately related to the e\ mixing parameter, is shown in Fig. ||. Here 
the present model is somewhat less repulsive than the various partial-wave analyses above 
150 MeV. 



The low-energy scattering parameters are shown in Table VIII and compared to exper 



imental results [p9| , p7tp2| . The scattering lengths and effective ranges are calculated both 



with and without the electromagnetic interaction. Without the electromagnetic interaction, 
the effective range function is simply given by F(k 2 ) = kcotS^ = — l/a+ ^rk 2 + 0(k 4 ). 
In the presence of the electromagnetic interaction, we have to use a more complicated ef- 



fective range function |p9|, where the phase shifts are with respect to the full long-range 



electromagnetic interaction. 

IV. PROJECTION INTO OPERATOR FORMAT 

We can project the strong interaction potential given above from S, T, T z states into an 
operator format with 18 terms 



E «pfa)0& • (25) 

p=l,18 



Here the first fourteen operators are the same charge-independent ones used in the Argonne 
v 14 potential and are given by 

Of" 1,14 = 1, Ti-Tj, (Ti-CT j,{(Ti-a j)(Ti-Tj), Sij,Sij(Ti-Tj), L • S, L- S(TyTj) , 

L 2 ,L 2 (Ti-Tj), L 2 ((7i-crj), L 2 ((7i-(Tj)(Ti-Tj), (L-S) 2 ,(L-S) 2 (Ti"0) • (26) 
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These fourteen components are denoted by the abbreviations c, r, a, or, t, tr, Is, 1st, 12, 
I2r, 12a, I2ar, ls2, and Is2t. The four additional operators break charge independence and 
are given by 

Ofj ' = Tij, (cr j • <t j)Tij , SijTij, (j~ z i + T z j) , (27) 

where Tij = 3r zi T z j — ti-t j, is the isotensor operator, defined analogous to the SV,- operator. 
These terms are abbreviated as T , aT, tT, and tz. The T, oT , and tT operators are charge- 
dependent and are "class II" forces, while the tz operator is charge-asymmetric and is a 
"class III" force (25|. 

The operator potential terms, v p , can be obtained from the channel potentials, Vg TNN , 
by a simple set of projections. We first introduce charge splitting for the central T = 1 
states, 

VS1,NN = 4l + V tl T H + V Tl(^i + Tzj) • (28) 

For the charge-independent potential this implies 

< = l Ml, PP + Vknn + ^Sl,np) ■ (29) 

We then project 



V c = 


U 9v ii 




(30a) 


V T = 




- 3< + - < ) , 


(30b) 


Va = 




+ vf - 3< - v«) , 


(30c) 


Var = 




- < - < + O , 


(30d) 



where of course v" = v^ np and Vq = v^ Q np . A similar set of projections is used for the L 2 
parts of the interaction. For the tensor, spin-orbit, and quadratic spin-orbit pieces, which 
exist only in S = 1 channels, the projections are (x = t, Is, ls2) 

^ = |0Ki + O , (31a) 
^ = |Ki-<o) • (31b) 
14 



The charge-dependent terms in Eq. (|2q) are given by 



V S1 — §[§( v Sl,pp + v Sl,nn) V Sl,np] > (32) 



51 — 6L2\ y Sl,H) T u Sl,nnJ u Sl,npl > 

which can be projected as 



t* = i(3wg + t;£) , (33a) 



«or = i«-<). (33b) 
The charge-dependent tensor term comes only from the spin-triplet channel, and reads 

V tT = \{Wn,pp + V U,nn) ~ v U,np] • ( 34 ) 

Finally, the charge-asymmetric terms are given by 

V S1 = \( v Sl,pp ~ v Sl,nn) ' (35) 

which leads to 

«™ = K 3u u + u £) ' ( 36a ) 

= J(«S- W S) • ( 36b ) 

As discussed in the previous section, we fix to reproduce the singlet nn scattering length 
by adjusting the parameter Pgi,nn ^° be slightly different from Pqi iPP - We are unaware of any 
nn data that would allow us to fix v^, but there have been numerous theoretical predictions 
for charge-symmetry breaking based on p-uo and n-n-n' mixing. Such models suggest that vH 
should be somewhat larger than VqI, but with a similar shape |33]. In the present work we 



make the simple assumption = v^f, which implies there is no v UTZ term. We also neglect 
the possibility of a charge-asymmetric tensor term v trz , which is why we end up with only 
one charge-asymmetric operator in our model. These choices are reflected in the parameters 
of Table |TJ. 

The first four operator components of the potential are shown in Fig. ^|. The tensor 
components are shown in Fig. [7| where we also show the CI part of the OPE potential used 
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here, and for comparison an OPE potential constructed using the same coupling constant 
and a dipole form factor (monopole at each nucleon-nucleon-pion vertex) with the cutoff 
mass A = 900 MeV. The spin-orbit and quadratic spin-orbit terms are shown in Fig. |8|, 
while the various L 2 components are shown in Fig. |9|. Finally, the charge-dependent and 
charge-asymmetric terms are shown in Fig. [It], along with the static Coulomb potential for 
comparison. 

The relative importance of the different CIB components is illustrated in Table [IX], where 
the evolution from the CI part of the interaction to the full pp interaction is displayed. The 
successive columns give the 1 Sq phase shifts for 1) the CI potential with an average nucleon 
and average pion mass, 2) with the correct proton mass, 3) with the correct CD OPE 
tail (i.e., correct neutral-pion mass) but the CI core, 4) with both the CD OPE and core 
interactions, and 5) with the electromagnetic potential added. From these it can be seen 
that the nucleon mass has a relatively small effect, while the CD OPE and core terms have 
relatively large effects at low energy, and the core contribution becomes dominant at higher 
energies. 



V. DEUTERON PROPERTIES 



The static deuteron properties are shown in Table [X] and compared to experimental 



values p4] , |3^ , ^|37| , p8| , |39l . The binding energy, E^, is fit exactly by construction. The 
expectation values for the kinetic energy, T, and for the EM, OPE, and remaining potentials 
are also shown. We note that the OPE potential dominates, while the EM potential gives 
a small but non-negligible 18 keV contribution, mostly from the magnetic moment term. 
The asymptotic S'-state normalization, As, the D/S-ratio, r], and the deuteron radius, 
rd, all come out close to the experimental values. The magnetic moment, /i^, and the 
quadrupole moment, Qj, are both underpredicted in impulse approximation; both have 
significant relativistic and meson-exchange corrections, as discussed below. Finally, the in- 
state percentage is about 5% smaller than that of the older Argonne Vu model [Ij] and almost 
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identical to that of the Paris potential 



The S- and D-wave components of the deuteron wave function are shown in Fig. [TT 
where they are compare to those for the older vu model. The short-range behavior of the 
wave function components is moderately different. The A(q 2 ) and B(q 2 ) structure functions 
and tensor polarization T 2 o(<2' 2 ) obtained with the present interaction model are displayed in 
Figs. lUll the experimental data is from Refs. |2|,|33iO>iIlil>ilSilil>ill • The model 



for the isoscalar electromagnetic current operator has been discussed in detail in Refs. |48| , f49 



here we only summarize its general structure, which consists of one- and two-body parts. 
The one-body part has the standard impulse approximation (IA) form, with inclusion, in the 
charge component, of the Darwin- Foldy and spin-orbit relativistic corrections |5(J. The two- 
body charge operators contain contributions that correspond (in an OBE picture) to those 
obtained from pion- and vector-meson (p and u) exchanges. These are obtained from the 
nonrelativistic reduction of the Born terms in the corresponding relativistic photoproduction 
amplitudes ||50| . The two-body current operators are constructed from the spin-orbit and 
quadratic momentum-dependent components of the interaction with the methods developed 
in Refs. We also consider the two-body charge and current operators associated 



with the p7T7 mechanism. In particular, we include in the nonrelativistic reduction of its 
current component the next to leading order correction arising from the tensor coupling 



of the p meson to the nucleon |49]. The Hohler parametrization 8.2 |)!| is used for the 
electromagnetic form factors of the nucleon, while an tu-pole term form factor is included at 
the p7T7 electromagnetic vertex. 

The calculated A(q 2 ) structure function is in excellent agreement with the experimental 
data over the whole range of measured momentum transfers. The Darwin-Foldy and spin- 
orbit relativistic corrections to the single-nucleon charge operator as well as the leading 
two-body charge contribution due to pion exchange play an important role, as it is evident 



from Fig. [12]. However, these same contributions lead to a significant discrepancy between 
theory and experiment in the tensor polarization. This observable and the A(q 2 ) structure 
function are mostly sensitive to the charge and quadrupole form factors. In particular, the 
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momentum transfer at which the minimum of T2o{<f) occurs is related to the position of the 
charge form factor zero. The relative shift between the predicted and experimental T 20 (g 2 ) 
minima implies, therefore, a corresponding shift between the charge form factor zeros. 

The calculated B(q 2 ) structure function is found to overpredict the experimental data in 
the momentum transfer range 10-45 fm -2 , and has a zero around 60 fm -2 . The leading two- 
body contributions are those due to the spin-orbit and quadratic spin-orbit components of 
the interaction. They are of opposite sign. However, the overestimate of the data indicates 
that the degree of cancellation between them is not quite enough. The pn^ current contri- 
bution is small over the momentum transfer range considered here (we have used the rather 
soft cutoff values of 0.75 GeV and 1.25 GeV at the irNN and pNN vertices, respectively, 



as suggested in Ref. [p3| .) At present, the two-body currents associated with the quadratic 
spin-orbit and L 2 components of the interaction are essentially obtained by minimal sub- 
stitution pi — > pi — [G s E (q 2 ) + Gg(g 2 )r 2j j]A(rj), where A is the vector potential, G% and 



G V E the isoscalar and isovector nucleon electric form factors J48|. It would be desirable to 
construct these current components in a more systematic way, as suggested in Ref. [Qj. 

Finally, the values for the quadrupole and magnetic moments obtained with the full 
charge and current operators are: Qd = 0.275 fm 2 and pd = 0.871 p Q . The measured 
quadrupole moment is underestimated by roughly 4%, while the measured magnetic moment 
is overestimated by 1.5%. The two-body charge and current contributions amount to 2% 
and 3% increases of the I A values for Q d and pd, as listed in Table |X|. 



VI. CONCLUSIONS AND OUTLOOK 

We have constructed a nonrelativistic NN potential with a local operator structure that 
gives an excellent fit to pp and np scattering data, as well as to low-energy nn scattering 
and the deuteron binding energy. We have projected the potential into charge- independent, 
charge-dependent, and charge- asymmetric pieces. In T = many-body systems only the 
CI part of the potential will contribute, while the CA part will contribute in systems with 
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T > | and the CD part for T > 1 systems. Because of the isotensor projection, the CI 
part automatically has the correct average of | pp (or nn) and \ np T = 1 interaction 
in the trinucleons, thus serving as a correct reference point for building NNN potentials. 
The CA part will contribute to the energy differences of mirror nuclei (the Nolen-Schiffer 
anomaly ||10|), while the CD part will contribute to the splitting of isobaric analog states, 
e.g., in the A = 6 nuclei. Studies of these effects are in progress. 

We have also computed the deuteron electromagnetic properties in both impulse ap- 
proximation and with relativistic and exchange- current corrections. The least satisfactory 
prediction of the potential is the small value for Qd, even after corrections are added. The 
full 4% discrepancy between the predicted and empirical Qd values is unlikely to be resolved 
by additional relativistic and/or two-body corrections not included in the present calcula- 



tion. Similar low values were found by the Nijmegen group in their recent fits [^4[ , which 
used rather different potential forms than the present model. We are less concerned about 
the 1.5% error in the magnetic moment or the overprediction of the B(q 2 ) structure function 
because of the uncertainties in the exchange currents discussed above. The A(q 2 ) structure 
function is very well reproduced, while the experimental tensor polarization T2o(q 2 ) still has 
rather large error bars above 10 fm~ 2 . 

Compared to the older Argonne vu potential, the present model has a weaker tensor 
force, which will generally lead to more binding in light nuclei and less rapid saturation in 
nuclear matter. This is counteracted by the weaker attraction in T = 1 NN states because 
of the mix of pp and np components. Initial calculations of few-body nuclei with the v is 
model show a slight net reduction in the binding energies of 3 H and 4 He compared to the 
v 14 model. Another feature of the new model is a moderately greater attraction in P waves. 
Few-body nuclei are not sensitive to this part of the AW interaction. However, preliminary 
calculations [fjS] of the binding energy of 16 O show a significant improvement in the relative 
stability of 16 and 4 He, which has been a persistent problem |TI|. We believe the Argonne 
Vig potential has an promising future for use in microscopic nuclear many-body theory. 
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TABLES 



TABLE I. Values of fundamental constants adopted in this work. 

he 197.32705 MeVfm 

m^o 134.9739 MeV/c 2 

7tv± 139.5675 MeV/c 2 

M p 938.27231 MeV/c 2 

M n 939.56563 MeV/c 2 

a- 1 137.03599 

Hp 2.79285 /i 

fi n -1.91304 hq 
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TABLE II. Short-range potential parameters in MeV. The asterisk denotes that the value was 



computed by Eq. (|23| ) and not fit. The three shape parameters are: c = 2.1 fm 2 , tq = 0.5 fm, 
and a = 0.2 fm. 



Channel 


Type 


I 


P 


Q 


R 


S = T = 1(W) 




-11 27028 

1 1 ■ | V / v. ,J 


3346 6874 


1859 5627* 


o 


S = 0,T = l(np) 


c 


-10.66788 


3126.5542 


1746.4298* 





S = 0,T = l(nn) 


c 


-11.27028 


3342.7664 


1857.4367* 





S = 0, T = 1 


Z2 


0.12472 


16.7780 


9 0972* 

■ y . vy *y i — i 


o 


S = 0, T = 




-2 09971 


1204.4301 


511 9380* 

•_y iiit/ ".y v. / v / i 


o 




/2 


-0.31452 


217 4559 

— .1. i • a '.y *y »y 


117 9063* 

j — l i • ■ y v / v / vy i 


o 


<? = 1 T = 1(W) 




-7.62701 


1815 4920 

J. l_ y 1 t. ) ■ x ».y V, 7 


969 3863* 


1847 8059 


5 = l,T = l(np) 


c 


-7.62701 


1813.5315 


966.2483* 


1847.8059 


5= 1,T = l(nn) 


c 


-7.62701 


1811.5710 


967.2603* 


1847.8059 


5 = 1,T = 1 


/2 


06709 


342 0669 


185 4713* 


-615 2339 




f 


1 07985 


o 


-190 0949 

i ■ y v / . v / ' y a '.y 


-811.2040 




Ls 


-0 62697 


-570 5571 


-309 3605* 


819 1222 

ky 1 • / . 1 6^ ^ ^ 




ls2 


0.74129 


9.3418 


5 0652* 


-376 4384 


S = 1,T = 


c 


-8.62770 


2605.2682 


1459.6345* 


441.9733 




12 


-0.13201 


253.4350 


137.4144* 


-1.0076 




t 


1.485601 





-1126.8359 


370.1324 




Is 


0.10180 


86.0658 


46.6655* 


-356.5175 




ls2 


0.07357 


-217.5791 


-117.9731* 


18.3935 
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TABLE III. Distribution of \ 2 by laboratory kinetic energy of the Nijmegen combined par- 



tial-wave analysis 


|| (PWA93) and the new Ar; 


^onne v is potential. Npp 


(N np ) denotes the number 


of pp (np) data in 


each enerj 


ry bin. 














x 2 {pp) 






x 2 


(np) 


Bin (MeV) 


Npp 


PWA93 




N np 


PWA93 


V18 


0.0-0.5 


134 


134.5 


136.3 


10 


9.7 


11.8 


0.5-2 


63 


39.7 


41.1 


5 


3.8 


7.4 


2-8 


48 


45.0 


36.0 


55 


52.4 


51.0 


8-17 


108 


103.0 


111.6 


182 


168.3 


164.8 


17-35 


59 


63.1 


72.2 


293 


226.6 


234.9 


35-75 


243 


213.4 


251.5 


328 


335.2 


339.3 


75-125 


167 


169.5 


171.5 


232 


237.1 


231.3 


125-183 


343 


379.7 


415.7 


333 


336.8 


363.5 


183-290 


239 


285.9 


304.8 


517 


494.6 


574.0 


290-350 


383 


360.7 


421.3 


571 


599.0 


708.0 


0-350 


1787 


1794.5 


1962.0 


2526 


2463.5 


2685.8 
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TABLE IV. pp phase shifts in degrees. Energies are in MeV. The 1 Sq includes the full electro- 
magnetic interaction (v EM (pp)) and is with respect to electromagnetic wave functions. The non-S 1 
waves are nuclear phase shifts of the Coulomb interaction including the form factor (Vci) with 
respect to Coulomb wave functions. 







'D 2 


3 Po 


3 Pi 


3 P 2 


£2 


3 F 2 


3 *3 


1 


32.68 


0.00 


0.14 


-0.08 


0.01 


-0.00 


0.00 


-0.00 


5 


54.74 


0.04 


1.61 


-0.90 


0.22 


-0.05 


0.00 


-0.01 


10 


55.09 


0.17 


3.80 


-2.05 


0.66 


-0.20 


0.01 


-0.03 


25 


48.51 


0.71 


8.78 


-4.89 


2.49 


-0.83 


0.10 


-0.23 


50 


38.78 


1.73 


11.75 


-8.23 


5.79 


-1.77 


0.32 


-0.69 


100 


25.01 


3.84 


9.61 


-13.11 


10.98 


-2.78 


0.73 


-1.47 


150 


15.00 


5.77 


4.72 


-17.27 


14.14 


-3.02 


1.06 


-1.96 


200 


6.99 


7.37 


-0.50 


-21.16 


15.91 


-2.88 


1.24 


-2.25 


250 


0.23 


8.61 


-5.50 


-24.86 


16.77 


-2.58 


1.21 


-2.45 


300 


-5.64 


9.52 


-10.17 


-28.37 


17.01 


-2.23 


0.90 


-2.66 


350 


-10.86 


10.14 


-14.49 


-31.70 


16.81 


-1.88 


0.29 


-2.95 
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TABLE V. nn phase shifts in degrees. Energies are in MeV. The 1 Sq includes the full electro- 
magnetic interaction {v EM (nn)). The non-S" waves only include the nuclear interaction. All phase 
shifts are with respect to Riccati-Bessel functions. 





'So 


X D 2 


3 Po 


3 Pi 


3 P 2 


82 


3 F 2 


3 F 3 


1 


57.07 


0.00 


0.21 


-0.12 


0.02 


-0.00 


0.00 


-0.00 


5 


60.64 


0.05 


1.88 


-1.04 


0.27 


-0.06 


0.00 


-0.01 


10 


57.48 


0.18 


4.17 


-2.24 


0.76 


-0.22 


0.01 


-0.04 


25 


48.80 


0.74 


9.13 


-5.12 


2.69 


-0.86 


0.11 


-0.24 


50 


38.47 


1.79 


11.89 


-8.48 


6.08 


-1.80 


0.32 


-0.70 


100 


24.45 


3.92 


9.48 


-13.38 


11.31 


-2.79 


0.74 


-1.49 


150 


14.38 


5.87 


4.46 


-17.58 


14.45 


-3.00 


1.07 


-1.98 


200 


6.34 


7.48 


-0.81 


-21.49 


16.19 


-2.84 


1.25 


-2.26 


250 


-0.42 


8.72 


-5.85 


-25.21 


17.00 


-2.53 


1.21 


-2.46 


300 


-6.31 


9.62 


-10.54 


-28.73 


17.20 


-2.18 


0.88 


-2.68 


350 


-11.53 


10.24 


-14.87 


-32.08 


16.96 


-1.83 


0.25 


-2.97 
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TABLE VI. np T = 1 phase shifts in degrees. Energies are in MeV. The 1 S includes the full 
electromagnetic interaction (v EM (np)). The non-5 waves only include the nuclear interaction. All 
phase shifts are with respect to Riccati-Bessel functions. 





'So 




3 Po 


3 Pi 


3 P2 


£2 


3 F 2 


3 F 3 


1 


62.02 


0.00 


0.18 


-0.11 


0.02 


-0.00 


0.00 


-0.00 


5 


63.50 


0.04 


1.64 


-0.93 


0.26 


-0.05 


0.00 


-0.00 


10 


59.78 


0.16 


3.71 


-2.04 


0.72 


-0.19 


0.01 


-0.03 


25 


50.61 


0.68 


8.32 


-4.82 


2.57 


-0.77 


0.08 


-0.20 


50 


40.09 


1.70 


10.99 


-8.15 


5.86 


-1.68 


0.28 


-0.61 


100 


26.02 


3.81 


8.69 


-13.07 


11.00 


-2.69 


0.67 


-1.35 


150 


15.98 


5.72 


3.78 


-17.28 


14.12 


-2.95 


0.98 


-1.82 


200 


8.00 


7.30 


-1.43 


-21.22 


15.86 


-2.82 


1.15 


-2.10 


250 


1.28 


8.52 


-6.41 


-24.95 


16.70 


-2.54 


1.10 


-2.30 


300 


-4.54 


9.43 


-11.06 


-28.49 


16.91 


-2.21 


0.77 


-2.51 


350 


-9.71 


10.06 


-15.36 


-31.85 


16.69 


-1.88 


0.14 


-2.81 
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TABLE VII. np T = phase shifts in degrees. Energies are in MeV. The coupled 3 S 1 - 3 D 1 
channel includes the full electromagnetic interaction (v EM (np)). The non-5 waves only include 
the nuclear interaction. All phase shifts are with respect to Riccati-Bessel functions. 





'Pi 


^3 


3 Si 


£\ 


3 D 1 


3 D 2 


3 D 3 


£3 


3 G 3 


1 


-0.19 


-0.00 


147.75 


0.11 


-0.00 


0.01 


0.00 


0.00 


-0.00 


5 


-1.51 


-0.01 


118.18 


0.66 


-0.17 


0.22 


0.00 


0.01 


-0.00 


10 


-3.11 


-0.07 


102.62 


1.14 


-0.65 


0.85 


0.01 


0.08 


-0.00 


25 


-6.48 


-0.42 


80.68 


1.77 


-2.72 


3.71 


0.08 


0.55 


-0.05 


50 


-9.85 


-1.13 


62.89 


2.11 


-6.28 


8.94 


0.40 


1.61 


-0.26 


100 


-14.20 


-2.22 


43.51 


2.52 


-12.04 


17.10 


1.61 


3.50 


-0.93 


150 


-17.68 


-2.98 


31.19 


2.96 


-16.39 


21.85 


2.92 


4.88 


-1.74 


200 


-20.79 


-3.61 


21.94 


3.43 


-19.82 


24.20 


4.00 


5.88 


-2.58 


250 


-23.65 


-4.22 


14.45 


3.92 


-22.59 


25.06 


4.76 


6.61 


-3.41 


300 


-26.28 


-4.87 


8.13 


4.43 


-24.83 


25.01 


5.21 


7.16 


-4.20 


350 


-28.71 


-5.59 


2.65 


4.95 


-26.65 


24.41 


5.39 


7.59 


-4.96 
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TABLE VIII. Scattering lengths and effective ranges in fm. 



Experiment 



Argonne v±s 



w/o V 



EM 



(Ipp 


-7.8063±0.0026 a 


-7.8064 


-17.164 


1 r 


2.794 


± 


0.014 a 


2.788 


2.865 


&nn 


-18.5 


± 


0.5 b 


-18.487 


-18.818 


T-nn 


2.8 


± 


0.1 b 


2.840 


2.834 


0"np 


-23.749 


± 


0.008° 


-23.732 


-23.084 


' np 


2.81 


± 


0.05 c 


2.697 


2.703 


3 a 

"-np 


5.424 


± 


0.003 c 


5.419 


5.402 


^np 


1.760 


± 


0.005 c 


1.753 


1.752 



a Ref. H| 
b Ref. |7| 
c Ref. |2] 
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TABLE IX. Evolution of x 5q pp phase shifts from the charge-independent potential to the full 



interaction, 


as described 


in the text. Energies are 


in MeV. 








CI 


+ m p 


+ CD v* 


+ CD v R 


+ V EM 


1 


58.40 


57.20 


57.42 


55.50 


32.68 


5 


61.44 


61.34 


60.88 


59.78 


54.74 


10 


58.14 


58.07 


57.71 


56.84 


55.09 


25 


49.33 


49.29 


49.05 


48.36 


48.51 


50 


38.95 


38.93 


38.76 


38.13 


38.78 


100 


24.93 


24.92 


24.80 


24.19 


25.01 


150 


14.87 


14.87 


14.77 


14.16 


15.00 


200 


6.86 


6.86 


6.77 


6.15 


6.99 


250 


0.11 


0.12 


0.04 


-0.60 


0.23 


300 


-5.75 


-5.74 


-5.82 


-6.47 


-5.64 


350 


-10.96 


-10.95 


-11.01 


-11.69 


-10.86 
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TABLE X. Static deuteron properties. 





Experiment 


Argonne v\% 


+ R + MEC 


Units 


E d 


2.224575(9) a 


2.224575 




MeV 


< T > 




19.814 




MeV 


^ EM - 
< V > 




0.018 




MeV 


<v 1T > 




-21.286 




MeV 


< v n > 




-0.770 




A If \T 

MeV 


A 


0.8846(8) 


0.8850 




£ 1/2 

fan ' 


rj 


0.0256(4) L 


0.0250 






I'd 


1.9660(68) d 


1.967 




fm 




0.857406(l) e 


0.847 


0.871 


IM) 


Qd 


0.2859(3) f 


0.270 


0.275 


fm 2 


Pd 




5.76 




% 



a Ref. @ 
b Ref. |5] 
c Ref. 
d Ref. |7] 
e Ref. |8] 
f Ref. M 
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FIGURES 

FIG. 1. Form factors in the electromagnetic interaction. 

FIG. 2. Phase shifts in the 1 S'o channel for np, nn, and pp scattering, compared to various 
partial-wave phase-shift analyses. 

FIG. 3. Phase shifts in the 3 Pq channel for np, nn, and pp scattering, compared to various 
partial-wave phase-shift analyses. 

FIG. 4. The e\ mixing parameter compared to various partial-wave phase-shift analyses. 

FIG. 5. Phase shifts in the l P\ channel, compared to various partial-wave phase-shift analyses. 

FIG. 6. Central, isospin, spin, and spin-isospin components of the potential. The central 
potential has a peak value of 2031 MeV at r = 0. 

FIG. 7. Tensor and tensor-isospin parts of the potential. Also shown are the OPE contribution 
to the tensor-isospin potential, and for comparison, an OPE potential with a monopole form factor 
containing a 900 MeV cutoff mass. 

FIG. 8. Spin-orbit and quadratic spin-orbit components of the potential. 

FIG. 9. I? components of the potential. 

FIG. 10. Charge-dependent and charge-asymmetric components of the potential. Also shown 
for comparison is the Coulomb potential, Vciipp)- 

FIG. 11. The deuteron S- and .D-wave function components divided by r. 
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FIG. 12. The deuteron electromagnetic structure function A(q 2 ) in impulse approximation 
(dashed line) and with relativistic and exchange-current corrections (solid line). Data are from 



Bonn @, Mainz |7j], Saclay g|, and SLAC g|. 



FIG. 13. The deuteron electromagnetic structure function B{q 2 ) in impulse approximation 
(dashed line) and with exchange-current corrections (solid line). Data are from Bonn [40], 



Mainz |3J, Saclay and SLAC |3|. 



FIG. 14. The deuteron tensor polarization T20 in impulse approximation (dashed line) and 
with relativistic and exchange-current corrections (solid line). Data are from Bates |4|],^5| and 



Novosibirsk |4| j47j] . 
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